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Abstract 


Infectious bronchitis coronavirus (IBV), Newcastle disease virus (NDV), and avian 
influenza virus (AIV) H9 subtype are major pathogens of chickens causing serious 
respiratory tract disease and heavy economic losses. To better understand the replication 
features of these viruses in their target organs and molecular pathogenesis of these different 
viruses, comparative proteomic analysis was performed to investigate the proteome changes 
of primary target organ during IBV, NDV and AIV H9 infections, using 2D-DIGE followed 
MALDI-TOF/TOF-MS. In total, 44, 39, 41, 48 and 38 proteins were identified in the tracheal 
tissues of the chickens inoculated with IBV (ck/CH/LDL/97I, H120), NDV (La Sota), and 
AIV H9, and between ck/CH/LDL/97]I and H120, respectively. Bioinformatics analysis 
showed that IBV, NDV and AIV H9 induced similar core host responses involved in 
biosynthetic, catabolic, metabolic, signal transduction, transport, cytoskeleton organization, 
macromolecular complex assembly, cell death, response to stress and immune system process. 
Comparative analysis of host response induced by different viruses indicated differences in 
protein expression changes induced by IBV, NDV, and AIV H9 may be responsible for the 
specific pathogenesis of these different viruses. Our result reveals specific host response to 
IBV, NDV, and AIVH9 infections and provides insights into the distinct pathogenic 


mechanisms of these avian respiratory viruses. 
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1 Introduction 


Avian respiratory viruses, including infectious bronchitis coronavirus (IBV), Newcastle 
disease virus (NDV), and avian influenza virus (AIV), are the primary cause of morbidity and 
mortality worldwide in the poultry industry and have been identified as the most 
economically important etiological agents of acute, highly contagious poultry diseases. 
Among these viruses, NDV is a member of the Avulavirus genus within the Paramyxoviridae 
family and is prevalent worldwide and infects almost all avian species. In fact, highly 
contagious respiratory, enteric, and neurological diseases in chickens are often caused by a 
variety of infectious NDV strains [1, 2]. AIV is an enveloped virus with a segmented, single- 
stranded, negative-sense RNA genome belonging to the Orthomyxoviridae family and can 
cause clinical diseases in a wide range of mammalian and avian species. Low pathogenic 
avian influenza virus subtype H9 circulates worldwide in multiple avian species thereby 
posing a potential threat to not only the poultry industry, but also to public health [3, 4]. IBV 
is a gammacoronavirus of the subfamily Coronavirinae, family Coronaviridae, and order 
Nidovirales [5], and primarily replicates in epithelial cells of the respiratory tract and causes 
lower-respiratory tract infections in chickens (Gallus domesticus). In addition, emerging IBV 
isolates that vary in their replicative abilities in non-respiratory tissues have been associated 
with proventriculus and kidney lesions [6, 7]. 

Proteomic analyses of pathogen-infected organisms can offer significant insight to the 
“battle” between pathogens and host during infection, hence proteomic analyses of host 


cellular responses against viral infections could aid in identification of potential host factors 
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involved in the viral life cycle and host defense mechanisms. To date, various proteomic 
strategies have been applied to elucidate protein dynamics during influenza infections, such 
as differential expression of brain proteins in response to the neurovirulent ATV H5N1 
infection [8], host responses of murine lung tissues infected with ATV H5N1 [9], and cellular 
protein changes in H1N1-infected human lung A549 cells [10]; while cellular responses to 
AIV H9N2 subtype infections has only been reported by Liu et al. in human cells [11]. For 
IBV, proteomic changes have been reported in the cytoplasm, nuclei, and nucleoli in 
mammalian Vero cells and nucleolar proteomes in avian cells (DF-1) infected with IBV 
(Beaudette strain) [12, 13]. We previously reported alterations in the proteome of tracheal and 
kidney cells during IBV strain H120 infection in vitro, and IBV strain ck/CH/LDL/97I and its 
attenuated strain ck/CH/LDL/97I P15 infection in vivo [14, 15]. However, alterations of 
cellular protein expression profiles, whether in NDV-infected cell lines or tissues, and the 
differences in induced protein expression changes between different avian respiratory viruses 
(IBV, NDV and AIV H9) have not yet been reported. 

The tracheal epithelium of chicken is not only the primary target for replication and 
infection of avian respiratory viruses, but also the first line of defense in the innate immune 
response. The avian respiratory viruses, including IBV, NDV and AIV H9, can infect tracheal 
epithelial cells and induce damage to tracheal epithelium resulting in respiratory disease. 
However, the severity and outcome of the clinical illness caused by these different viruses are 
different [16-18]. To better understand the molecular basis of pathogenesis and specific virus- 
host interactions, herein, we described the alterations of protein expression profiles in chicken 


tracheal tissues caused by IBV ck/CH/LDL/97I and H120, NDV La Sota, and ATV H9 
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infection by using 2D-DIGE analyses. A series of differentially expressed proteins associated 
with viral life cycle and host response to IBV, NDV and AIV H9 infections were identified. 
The present proteome results using the same virus of different strains (IBV virulent 
ck/CH/LDL/97I and H120 vaccine) and different viruses IBV, NDV and AIV H9) might 
provide more significant information on specific aspects of molecular pathogenesis and virus- 
host interactions of these different viruses. Furthermore, these findings might also provide 
useful clues for the development of novel prevention or therapeutics strategies against these 


viruses. 
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2 Materials and methods 


2.1 Viral strains 


A commercial vaccine against IBV strain H120 was obtained in freeze-dried 1000-dose 
vials that each contained at least 10° egg infective dose 50% (EIDs0) per dose (isbi® BIO 
H120; batch No. F48355). IBV strain ck/CH/LDL/97I [19], NDV strain La Sota vaccine, and 
AIV subtype H9 [20] were used in this study. The virus stock for this study was produced by 
inoculating the virus into embryonated SPF chicken eggs via the allantoic cavity and 
collecting the infectious allantoic fluid 72 h post-inoculation. The allantoic fluid was clarified 
by centrifugation at 3000 x g for 10 min and filtered with a Teflon membrane and the 
presence of IBV was confirmed by negative contrast electron microscopy (JEM-1200, EX; 
JEOL Ltd., Tokyo, Japan) and reverse transcription polymerase chain reaction (RI-PCR) as 
previously described [21]. NDV and subtype H9 AIV viral particles in the allantoic fluids of 
inoculated eggs were confirmed using hemagglutination inhibition (HD) analysis with specific 
antibodies, respectively [22]. The titers of the viruses were determined by inoculation at 10- 
fold dilutions into groups of five 10-day-old embryonated chicken eggs. The EID59 was 


calculated by the method of Reed and Muench [23]. 
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2.2 Experimental design 


One hundred and twenty 1-day-old SPF White Leghorn chicks were divided into five 
groups of 24 birds each and housed in different isolators with negative pressure. Chickens 
assigned to groups 1-4 were inoculated with IBV ck/CH/LDL/971, H120, NDV La Sota, and 
AIV H9, respectively, by oculonasal application at 3 days of age with a dose of log10° 
EIDs50/100 uL per chick. Birds in group 5 were mock-inoculated with sterile allantoic fluid as 
a negative control. This study was approved by the Animal Welfare Committee of 
Heilongjiang Province, China. All experiments were conducted at appropriate biosafety levels. 

Three birds from each group were killed humanely at 3, 7, 12, and 21 days post- 
inoculation (dpi) and fresh tissue of trachea samples were quickly harvested and washed with 
ice-cold Tris-sucrose buffer (10 mM Tris, 250 mM sucrose, pH 7.5), snap-frozen in liquid 
nitrogen, and stored at —80°C for subsequent 2D-DIGE analysis. Sera samples from 10 
chickens in each group were collected at 3, 7, 12, and 21 dpi for enzyme-linked 
immunosorbent assay (ELISA) and HI analyses. Sera samples from groups 1, 2, and 5 were 
assayed using a commercial IBV antibody test kit (IDEXX Corporation, Westbrook, Maine, 
USA) according to the manufacturer’s instructions. Serum samples from groups 3, 4, and 5 
were detected for HI antibody titer (positive when HI >3) against NDV and AIV, respectively, 
using an HI assay as reported previously [24]. Viral load of IBV, NDV and AIV in tracheal 
tissues were analyzed using methods described previously. Briefly, Viral RNA was extracted 
from the tracheal tissues of sacrificial birds from each group using TRIzol Reagent 


(Invitrogen, Carlsbad, CA, USA). For detection of the viruses, specific primers and TaqMan 
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probes against IBV [25], NDV [26], and AIV subtype H9 [4] viral RNA were selected. Real- 
time RI-PCR was performed on a LightCycler® 480 real-time PCR system with One Step 
PrimeScript™ RT-PCR Kit (TaKaRa, Dalian, China) [27]. All of the samples were tested in 
triplicate in each reaction. Standard templates with known concentration and no template 


negative control were contained for quantitative analysis. 


2.3 2D-DIGE analysis of tracheal proteins 


The tracheal tissues were ground into powder in liquid nitrogen. Fifty milligrams of 
powder were extracted with 500 uL of DIGE-compatible lysis buffer containing 7 M urea, 2 
M thiourea, 4% CHAPS, 1% nuclease mix, and 1% Protease Inhibitor Mix (GE Healthcare, 
Waukesha, WI, USA), vortexed for 3 h in an ice bath, and then centrifuged at 15 000xg for 1 
h at 4°C. The supernatant was collected and purified using the Ettan 2-D Clean-Up Kit (GE 
Healthcare). The pH of the protein sample was adjusted to 8.5 and the protein concentration 
was determined using the Ettan 2-D Quant Kit (GE Healthcare). The purified protein lysates 
(50 ug) were labeled with 400 pmol of Cy2, Cy3, and Cy5 dyes (GE Healthcare) following 
the manufacturer’s instructions. Briefly, the Cy2-labeled internal standard was created by 
pooling an aliquot from all involved biological samples, while Cy3 and Cy5 were used 
randomly to label samples from the IBV ck/CH/LDL/97]-, H120-, NDV La Sota-, AIV H9-, 
and mock-infected groups. The experimental design of the different fluorescent dye labeling 
protocols is listed in Supporting Information Table S1. The labeling reaction was performed 
as described in our previous study [15] and labeled samples were subsequently loaded onto 
24-cm, linear pH 3-10 IPG strips (GE Healthcare) for integrated isoelectric focusing (IEF) 
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which was performed using a Ettan™ IPGphor 3 fully integrated IEF system (GE Healthcare). 
After IEF, the IPG strips were treated in an equilibration buffer (6 M urea, 50 mM Tris-HCl, 
pH 8.8, 30% glycerol, 2% SDS, and 0.002% bromophenol blue) containing 1% dithiothreitol 
and 2.5% iodoacetamide. A second dimension sodium dodecyl sulfate polyacrylamide gel 
electrophoresis (SDS-PAGE) was performed using an Ettan DALT six system electrophoresis 
system (GE Healthcare). 

After SDS-PAGE, the CyDye-labeled gels were scanned using the Typhoon 9400 
scanner (GE Healthcare) at a resolution of 100 um. The intensity was adjusted to ensure that 
the maximum volume of each image was within 60,000—80,000. DIGE analysis was 
performed using DeCyder software version V6.5 (GE Healthcare). Protein spots were 
pairwise compared among control and different infection groups and differentially expressed 
spots with statistically significant results via the Student’s t-test (p < 0.05) and average ratio 
>1.2 or <-1.2 were selected for MS identification. Spot picking was performed using 
preparative gels under the conditions as described above for 2D-DIGE except that the IPG 
strips were loaded with 1000 ug of protein and the gels were stained with Coomassie Blue R- 
350, then the protein spots were excised manually from the gels and subjected to MALDI- 


TOF/TOF MS identification. 
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2.4 MALDI-TOF/TOF MS and database search 


The gel spots were excised, trypsin-digested, and subjected to MALDI-TOF/TOF MS 
analysis as described in our previous study [15]. The peptide mass fingerprint (PMF) 
combined MS/MS data were submitted to MASCOT version 3.0 (Matrix Science, Inc., 
Boston, MA, USA) for identification according to the National Center for Biotechnology 
Information (NCBI) nr database (release 16/01/2010, 10,343,571 sequences, and 
3,528,215,794 residues). The search parameters were set as follows: Gallus, trypsin cleavage 
(one missed cleavage allowed), carbamidomethylation of cysteine as a fixed modification, 
oxidation of methionine as a variable modification, peptide mass tolerance set at 100 ppm, 
and fragment tolerance set at 0.8 Da. The criterion for successful identification of a protein 


was a protein score confidence interval (CI %) >95%. 


2.5 Gene ontology (GO) annotation and classification 


To investigate the potential roles played by the differentially expressed proteins during 
infection, GI numbers of identified proteins were submitted to GORetriever 


(http://www.agbase.msstate.edu/) to obtain the Gene Ontology (GO) annotations. If no 


annotation was acquired from GORetriever, GOanna was used to retrieve GO annotations 
assigned on the basis of searches using the Basic Local Alignment Search Tool. Then, the 


resulting annotations were summarized using GOSlimViewer. 
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2.6 Analysis of mRNA abundance by real-time RT-PCR 


Ten target proteins which extensively changed during IBV, NDV and AIV H9 infections 
shown by 2D-DIGE and potentially involved in the viral infection were selected for 
transcriptional analysis. The expression of selected proteins were not all altered at four time 
points (Supporting Information Table S2), therefore the mRNA levels of selected proteins at 
their corresponding time point were analyzed, respectively. Specific primers for the target 
genes of MS/MS-identified proteins were designed according to the published gene 
sequences using Beacon Designer software 7.5 (Premier Biosoft International, Palo Alto, CA, 
USA) and information on the primers is listed in Supporting Information Table S3. Total 
tracheal tissue RNA was extracted from sacrificed birds of each group at each time point 
using TRIzol Reagent (Invitrogen) and then reverse-transcribed to cDNA. The real-time PCR 
was performed using a LightCycler® 480 real-time PCR system (Roche Diagnostics 
Deutschland GmbH, Mannheim, Germany). Melting curve and quantitative analyses of the 
data was determined using the LightCycler® 480 real-time PCR system with the relative 
quantification 2A" model, Cp values of target gene in control group at each time point were 
used as the calibrator (relative expression=1), and 28S ribosomal RNA gene was used as an 


internal reference gene. 
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2.7 Western blot analysis 


To verify the 2D-DIGE results, several representative proteins were selected for Western 
blot analysis based on the commercial availability of the corresponding antibodies. Tracheal 
proteins from sacrificed birds of each group at all four time points were extracted and 
equivalent lysate amounts were subjected to 12% SDS-PAGE and then transferred to 
nitrocellulose membranes. After blocking overnight at 4°C, the membranes were incubated 
with goat polyclonal antibody to annexin A2 (ANXA2; sc-1924; Santa Cruz Biotechnology), 
rabbit polyclonal antibody to annexin Al (ANXA1; sc-11387; Santa Cruz Biotechnology), 
goat polyclonal antibody to peroxiredoxin I (PRDX1; sc-23970; Santa Cruz Biotechnology), 
mouse monoclonal antibody to heat shock protein 27 (HSPB1; sc-51956, Santa Cruz 
Biotechnology), rabbit polyclonal antibody to phospho-HSPB1 (pSer82) (SAB4504434, 
Sigma-Aldrich) and mouse monoclonal antibody to B-tubulin (sc-55529; Santa Cruz 
Biotechnology) for 1 h at 37°C, respectively. Then, the membranes were incubated with 
IRDye800DX “conjugated anti-mouse, -goat, and -rabbit secondary antibodies (610-130-121, 
605-430-002, and 611-730-127; Rockland Immunochemicals, Gilbertsville, PA, USA), 
respectively, for 1 h at 37°C. After washing, the membranes were scanned using a LI-COR 


infrared imaging system with Odyssey software (Li-Cor Biosciences, Lincoln, NE, USA). 
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3 Results 


3.1 Clinical observations, serology and virus detection in trachea 

All chickens infected with the IBV ck/CH/LDL/97I strain exhibited clinical respiratory 
signs at ~3—12 dpi and one of the diseased birds died at 12 dpi. The clinical signs included 
tracheal rales, watery eyes, nasal mucus, and sneezing. Gross lesions within the chickens 
were confined mainly to the kidneys [15], whereas the chickens in the control group and 
those inoculated with the IBV H120 and NDV La Sota strains showed no apparent respiratory 
clinical signs and no gross lesions. ATV H9-infected chickens appeared lethargy, depression, 
coughing, sneezing and dyspnea during the experimental period. 

The serological responses induced by infections with IBV ck/CH/LDL/97I, H120, NDV 
La Sota, and AIV H9 are summarized in Supporting Information Table $4. No seroconversion 
was observed in any group at 3 dpi. Antibodies appeared at 7 dpi and all chickens showed 
seroconversion after 21 dpi in the IBV, NDV, and AIV infection groups. The chickens in the 
control group showed negative serum antibody responses against IBV, NDV, and AIV. 

IBV, NDV, and AIV infections were verified using real-time RT-PCR. As shown in Fig. 1, 
negative results were obtained following viral detection in the tracheal tissues from chickens 
in the control group. In the IBV ck/CH/LDL/97]-infected group, the IBV genome was 
detected at 3, 7, 12, and 21 dpi. The viral load gradually decreased during the experimental 
period after the viral copy number peaked at 3 dpi (Fig. 1A). Similar results were obtained 
from group H120, except no viral genome was detected at 21 dpi (Fig. 1B). In the NDV La 
Sota- and AIV H9-infected groups, viruses were detected at 3 dpi and viral copy number 
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peaked at 7 dpi, whereas the results at 12 and 21 dpi were negative (Fig. 1C and D). The 
differences in viral kinetics among IBV ck/CH/LDL/971, H120, NDV La Sota and AIV H9 
were demonstrated in Fig. 1E. Taken together, these results confirmed the successful virus 


infections to SPF chicken. 


3.2 2D-DIGE profiles of tracheal tissues following infections with avian respiratory 


viruses (IBV, NDV, and AIV H9) 


To obtain a detailed comparison of the differential protein expression profiles, the 
tracheal proteins from avian respiratory virus (IBV ck/CH/LDL/97I, H120, NDV La Sota, 
and AIV H9)-infected and mock-inoculated chickens were extracted for quantitative and 
comparative analysis via 2D-DIGE at 3, 7, 12, and 21 dpi. Comparative analysis revealed that 
total 53, 44, 70, and 33 spots were significantly differentially expressed at 3, 7, 12, and 21 dpi 
respectively, which are indicated by arrows in the corresponding representative 2D-DIGE 
map (Fig.2A, B, C and D). In the IBV ck/CH/LDL/97]-, H120-, NDV La Sota-, and ATV H9- 
infected groups, and between ck/CH/LDL/97I- and H120-infected groups, 72, 65, 63, 86 and 
63 differentially expressed protein spots were identified, respectively (Table 1 and Supporting 
Information Figure $1). In general, these data showed that the expression of more host 
protein spots were changed in the chicken tracheal tissues after virulent IBV ck/CH/LDL/97I 
and AIV H9 infections, comparing with those of the vaccine strains IBV H120 and NDV La 


Sota strains. 
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3.3 Identification and classification of differentially expressed proteins 


To identify the differentially expressed proteins in viral infected tracheal tissues, 200 
spots with thresholds >1.2-fold (Supporting Information Table S2) were excised from 
preparative gels and subjected to MALDI-TOF/TOF MS identification. Detailed information 
of the identified proteins is provided in Supporting Information Table S5 and Supporting 
Information Fig. S2. Totally, 200 differentially expressed protein spots corresponding to 75 
altered cellular proteins were successfully identified; the biological functions of the 
differentially expressed proteins were summarized in Supporting Information Table S6. 
Number of differentially expressed spots and proteins induced by IBV ck/CH/LDL/97I, H120, 
NDV La Sota, and AIV H9 infections were summarized in Table 1. Among these identified 
proteins, some were identified from two or more separated spots, indicating that they may be 
differently modified (e.g. phosphorylation of ANXA1 and HSPB1) in the process of host 
response to IBV, NDV and AIV H9 infections (Fig. 2). In order to verify the expression of 
phosphorylated HSPB1, tracheal proteins from sacrificed birds of each group were examined 
by Western blot analysis with specific antibodies to phospho-HSPB1 (pSer82). The up- 
regulation of phosphorylated HSPB1 in IBV-, NDV- and AIV H9-infected groups at 7dpi and 
12dpi were indentified (Fig. 3A and 3B); these results were consistent with the expression 
pattern of HSPB1 spots showed by 2D-DIGE, confirming the expression of phosphorylated 
HSPB1 during IBV, NDV and AIV H9 infections. 

To better understand the potential role of host factors involved in viral infection, 


differentially expressed proteins were classified based on gene ontology annotations derived 
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from the AgBase database (www.agbase.msstate.edu/). As illustrated in Fig. 4A, the biological 


process of proteins differentially expressed in IBV ck/CH/LDL/97I-, H120-, NDV La Sota-, 
and AIV H9- infected groups and between IBV ck/CH/LDL/97I- and H120-infected groups 
(Fig. 4B) all mainly involved in anatomical structure development, biosynthetic process, 
transport, signal transduction, protein complex assembly, cell death, cell differentiation, 
translation, cellular protein modification process, response to stress, cytoskeleton 
organization, growth, macromolecular complex assembly, immune system process, vesicle- 
mediated transport, embryo development, cell proliferation, transmembrane transport and 
various metabolic associated process. These analyses indicated that all of the four viruses 
affect host global response similarly. 

In order to highlight the differences of proteome changes in trachea induced by different 
viruses, a side-by-side comparison of the differentially expressed proteins in trachea between 
different groups was performed (Table 2 and Supporting Information Table S7). The 
comparisons revealed that alteration of certain proteins and corresponding biological process 
appeared to be general or specific to host response in the process of different virus infections. 
For example, ENO1, ENO3, PGAM1, ANXA1, CKB, GST, HSPB1, HBAA, HBAD, LMNA, 
akr and PRDX1 were commonly altered in all viruses incubated groups, likely suggesting that 
these proteins were involved in and had similar functions in response to the infection of the 
three avian respiratory viruses. GAPDH, LDHA, PPIA, PRDX6 and HADH were altered in 
pathogenic IBV ck/CH/LDL/97I- and AIV H9- infected groups, while PSMB2 was more 
abundant in vaccine strains IBV H120- and NDV La Sota- incubated groups. And we also 


found that some proteins were specifically affected by IBV ck/CH/LDL/97I (S100A11, 
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ANXA8 and PHB), H120 (UQCREFS1), NDV La Sota (rbf and VDAC1) and AIV H9 (RPSA, 
CTHRC1, PSMA6 and TAGLN) infections, respectively (Table 2). 

Expectedly, some proteins identified in this study have also been reported in proteomic 
analyses of cells or tissues infected with IBV or AIV (Table 3). For example, up-regulation of 
LGALS1, PPIA, PRDX4 and HSPB | has been found in AIV H9- infected human cells [11]. 
ANXA2, LGALS1, EIF5A2 and PPIA were also up-regulated in IBV-infected avian DF-1 
cells [12]. And down-regulation of ANXA1 was reported in IB V-infected chicken trachea in 
our previous works [14, 15]. The common alteration of these proteins in different models of 
infection with the same virus of different strains suggests that these proteins certainly 


involves in the course of corresponding virus infection. 


3.4 Transcriptional profiles of differentially expressed proteins 


Transcriptional alteration analysis of ten selected genes at 3, 7, 12, and 21 dpi were 
achieved by quantifying the mRNA transcripts, in which the housekeeping gene 28S rRNA 
was used as an internal control (Fig. 5A). In this study, primers designed for transcriptional 
alteration analysis were located at positions across intron-exon boundaries to avoid the 
interference of genomic DNA, so quantitative comparison of total amount of mRNA 
transcripts of each gene was conducted. In general, the trends in the change in mRNA 
abundance of selected genes at four time points (Fig. 5A) were similar to the change patterns 
of their corresponding proteins showed by 2D-DIGE (Fig. 5B and Supporting Information 
Table S2). And degrees of change in mRNA abundance of ANXA1, HSPB1, SOD1, EIF5A2, 
LGALS1 and PRDX1 were more obvious than the abundance change of corresponding 
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proteins showed by 2D-DIGE. This is not unexpected because post-transcriptional 
mechanisms, including protein translation, post-translational modification, and degradation 
may also contribute to the abundance change of these proteins during virus infection [28]. 
These data provided transcriptional information complementary to the differentially 


expressed proteins revealed by proteomics analysis. 


3.5 Validation of the differentially expressed proteins by western blot analysis 


To confirm the dynamic alterations of protein abundance during ck/CH/LDL/97I, H120, 
La Sota, and H9 infections, four protein candidates ANXA1, PRDX1, ANXA2 and HSPB1 
were chosen for western blot analysis with specific antibodies (Fig. 3A and 3B). The 
abundance changes of ANXA1 and PRDX1 (both altered at all four time points); ANXA2 
(altered only at 12dpi) and HSPB1 (altered at 3,77, and 12dpi) were consistent with the 
expression changes showed by 2D-DIGE (Fig. 5B and Supporting Information Table S2). 
Altogether, Western blot validated the results obtained from proteomics analyses during IBV, 


NDV and AIV H9 infections. 
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4 Discussion 


In this study, host responses against different avian respiratory viruses which result in 
distinct disease phenotypes were analyzed using comparative proteomics. IBV 
ck/CH/LDL/97I and AIV H9 are two kinds of pathogenic viruses for chicken causing 
significantly clinical signs, while, IBV H120 and NDV La Sota are avirulent vaccine strains 
inducing effective immune response. Our proteomic analyses of the alterations in protein 
expression induced by different viruses under the same physiological conditions providing 
more specific information of molecular pathogenesis and virus-host interactions of these 
different viruses. Although GO functional annotation and classification showed that all four 
viruses affect host global cellular processes similarly, the trends and degrees of alteration of 
protein expression were not identical. What apparent from the results is that pathogenic IBV 
ck/CH/LDL/97I and AIV H9 induced slightly more changes in protein expression than the 
vaccine strains IBV H120 and NDV La Sota in chicken trachea. In addition, the differences in 
protein expression changes induced by IBV ck/CH/LDL/97]I, H120, NDV La Sota, and AIV 
H9 indicated that different interactions likely occurred in the host response to infection by 
these viruses and resulted in the different pathogenicity. This is the first report of comparative 
proteomics analysis performed to study the host response to infections of different avian 
respiratory viruses; the information obtained herein might provide critical clues for future 


studies. 
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4.1 Inflammatory responses in IBV, NDV and AIV H9 infections 


IBV, NDV and AIV H9 subtype viruses can initially infect the chicken tracheal 
epithelium and elicit inflammatory responses associated with the disease [16, 29, 30]. In this 
study, several proteins involved in inflammatory responses were identified including ANXA1, 


PPIA, LGALS1 and PRDX4. 


Our previous works have been reported down-regulation of the annexin Al (ANXA1) in 
IBV-infected trachea [14, 15]. In this study, the endogenous anti-inflammatory protein 
ANXAI [31] was found down-regulated in both pathogenic IBV ck/CH/LDL/97]- and AIV 
H9-infected groups at 7dpi and in all virus-infected groups at 12dpi, suggesting that the 
decreased expression of ANXA1 may be associated with the inflammatory phenotype during 
the course of avian respiratory viruses infections. More notably, ANXA1 was identified in 
more than one spots in our study, suggesting that ANXA1 phosphorylation during IBV, NDV 
and AIV H9 infections may associates with their pathogenesis [32]. The peptidyl-prolyl 
isomerase (PPIA) has been shown to elicit inflammatory responses and inhibit CD4* T cell 
activation through its interaction with Itk [33]. We noticed that pro-inflammatory PPIA was 
only up-regulated in pathogenic IBV CK/CH/LDL/97I- and AIV H9-infected groups; there 
was no alteration of PPIA expression in trachea of chickens infected with IBV H120 and 
NDV La Sota. This difference may indicate that pathogenic IBV ck/CH/LDL/97I and AIV H9 
trigger more severe inflammatory responses and immune suppression than IBV H120 and 
NDV La Sota vaccine strains, which associated with their different pathogenetic processes. In 
addition, PPIA can facilitate replication of human immunodeficiency virus-1 (HIV-1) [34] 
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and coronaviruses including IBV [35]. However, PPIA can inhibit the influenza A viral 
replication through enhancing the degradation of M1 protein [36]. Chicken PPIA also showed 
anti-influenza virus activity [37]. Up-regulated PPIA has been detected in human cells 
infected with AIV H9N2 [11]. Hence, we speculate that the up-regulation of PPIA may 
contribute to the replication and pathogenesis of IBV ck/CH/LDL/97]I and participate in host 
defense against AIV H9. 

Peroxiredoxin 4 (PRDX 4) can protect pancreatic beta-cells against injury by 
suppressing oxidative stress and inflammatory signaling activation [38]. Galectin-1 (LGALS1) 
is an important contributor to immune homeostasis and exhibits anti-inflammatory activity 
during infection [39]. We found that anti-inflammatory PRDX4 and LGALS1 were generally 
up-regulated in IBV ck/CH/LDL/97I-, NDV La Sota- and AIV H9-infected groups at 12dpi. 
In the case of avian respiratory viral infections, the inflammatory process diminished at the 
recovery stage [16]. Hence, the up-regulation of anti-inflammatory PRDX 4 and LGALS1 
may be involved in tracheal epithelium recovery from inflammation damage induced by IBV 
ck/CH/LDL/97I1, NDV La Sota and AIV H9. Besides as an immunomodulator, galectin-1 can 
function as pattern recognition receptor (PRR) to facilitate pathogen recognition [39]. The 
interaction of galectin-1 with viral envelope glycoproteins can inhibit infectivity of ATV and 
Nipah virus [40, 41]. Galectin-1 promotes HIV-1 infection through stabilization of viral 
adsorption [42]. Moreover, up-regulation of galectin-1 was also found in human cells infected 
with AIV H9N2 [11], suggesting that galectin-1 may participate in innate immune response to 
AIV H9 infection. Our results also indicate that galectin-1 potentially involved in IBV 


ck/CH/LDL/97I and NDV La Sota infections, whether galectin-1 binds to the surface of IBV 
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and NDV, and functions as effective PRR or promotes infection by facilitating viral 


adsorption need to be further elucidated. 


4.2 IBV, NDV and AIV H9 may induce alterations of nuclear structure and hijack 


host translational machinery 


Nucleus associated proteins, including nuclear mitotic apparatus protein (NUMA), 
lamin-A (LMNA), and eukaryotic initiation factor SA-1,-2 (eIF5A1, 2), were identified in the 
present study. Though the primary site of replication of IBV, NDV and AIV is the cytoplasm, 
IBV nucleocapsid protein, NDV matrix protein, and AIV nucleoprotein can localize to the 
nucleolus [43-45]. In this study, LMNA was universally down-regulated in all virus-infected 
groups and NUMA was down-regulated in IBV ck/CH/LDL/971, H120- and AIV H9-infected 
groups. NUMA and LMNA are components of nuclear matrix and nuclear lamina, down- 
regulations of NUMA and LMNA might implicate that IBV, NDV and AIV H9 infections 
commonly result in disruption of nuclear architecture due to nuclear entry and accumulation 
of viral proteins in the nucleolus. Down-regulation of LMNA has also been detected in 
influenza virus-infected cells [10] and IBV-infected trachea [15]. 

As well known that virus may hijack host translational machinery by regulating distinct 
elFs [46-49]. Previous reports showed eIFSA can regulate translation and replication of HIV- 
1 [48], and eIF5A mutants inhibit HIV-1 replication [49]. In this study, eIF5A1 was down- 
regulated in IBV H120-, NDV La Sota- and AIV H9-infected groups, while, eIF5A2 was up- 
regulated in IBV ck/CH/LDL/97I-, NDV La Sota- and AIV H9-infected groups. We speculate 
that IBV, NDV and AIV H9 may interfere with the host protein synthesis by down-regulating 


This article is protected by copyright. All rights reserved 


www.proteomics-journal.com Page 25 Proteomics 


constitutively expressed eIF5A1. While the up-regulation of eIF5A2 facilitates viral mRNAs 
translation in infected cells; the expression of eIF5A2 is extremely low in normal conditions. 
Furthermore, the differences in expression of eIF5A1 and eIF5A2 indicate that the two IBV 
strains differentially hijack host translational apparatus for their own benefits. Our results 
provide additional evidences that RNA viruses can interact with the nucleolus to hijack 


cellular functions and recruit nucleolus proteins to facilitate virus replication. 


4.3 IBV, NDV and AIV H9 manipulate host UPP through different strategies 


The ubiquitin-proteasome pathway (UPP), a cellular ubiquitin-dependent protein 
degradation system, has been implicated in viral infection cycle and virus-host interactions 
[50]. Inhibition of the UPP has been reported to affect the replication of several viruses, 
including influenza virus [51] and coronavirus [52]. Moreover, the interaction of HIV-1 Tat 
protein with proteasomal f subunits contributes to the immune escape of HIV-1 [53]. 
However, the potential roles of host UPP during IBV and NDV infections and precise 
interaction mechanisms between UPP and AIV are not well elucidated. Previous reports 
showed that proteasome subunit alpha and beta types were up-regulated in IBV-infected DF-1 
cells [12]. In this study, proteasomal subunit alpha type-6 (PSMA6) was up-regulated only in 
AIV H9-infected group, whereas proteasomal subunit beta type-2 (PSMB2) was commonly 
up-regulated in IBV H120- and NDV La Sota-infected groups and more abundant in H120- 
infected group than ck/CH/LDL/97]I-infected group, indicating that different viruses (IBV, 
NDV, and AIV) and different strains of the same virus (ck/CH/LDL/97I and H120) may 
evolve different strategies to utilize UPP to facilitate replication and infection by specific 
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interaction with distinct proteasomal subunits. 


4.4 Different roles of ANXA2 in IBV, NDV and ATV H9 infections 


Annexin A2 (ANXA2) has been shown involved in exocytosis, cell proliferation, 
membrane trafficking and apoptosis [54], the organized ANXA2 membrane domains are 
necessary for virion budding from infected cells [55]. Notably, many viruses utilize ANXA2 
for effective infection and replication, such as respiratory syncytial virus (RSV) [56], HIV 
[57] and HCV [58]. In this study, ANXA2 was commonly up-regulated during IBV H120, 
NDV La Sota and AIV H9 infections. We speculate that the up-regulated ANXA2 may 
support replication of AIV H9 by incorporated into influenza virus particles to enhance 
influenza virus replication via mediated activation of plasminogen [59]. However, up- 
regulated ANXA2 might act as an anti-viral protein to inhibit replication of IBV H120 as it 
was reported that ANXA2 reduces viral frameshifting efficiency and successful replication of 
IBV as a cellular RNA-binding protein [54]. Up-regulation of ANXA2 has also been found in 
avirulent IB V-infected DF-1 cells [12]. In our previous work, the abundance of ANXA2 was 
only increased following attenuated IBV ck/CH/LDL/97I P15 infection, no significantly 
abundance changes of ANXA2 were observed during virulent IBV ck/CH/LDL/97I Ps; 
infection, suggesting that the ANXA2 may be specific host response to avirulent IBV strains 
(ck/CH/LDL/97I Ps and H120). The function of ANXA2 in life cycle of NDV needs further 


investigation. 
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4.5 IBV, NDV and AIV H¢9 infections and stress response proteins 


Heat shock protein 27 (HSPB1) is a multifunctional protein involves in various cellular 
processes and viral infection [60]. Phosphorylated HSPB1 has been shown to participate in 
the replication of adenovirus [60], HCV [61], and herpes simplex virus type 1 [62]. In this 
study, the up-regulation of phosphorylated HSPB1 during IBV, NDV and AIV H9 infections 
was confirmed by western bolt (Fig. 3). It was reported that RSV induced phosphorylation of 
HSPB1 in human bronchial epithelial cells ultimately results in epithelial membrane 
permeability [63]. Avian respiratory viruses infections can result in tracheal epithelium lesion, 
but the severity and persistence of the lesions differed among different viruses [16-18]. 
Therefore, phosphorylated HSPB1 which mainly up-regulated in IBV H120- and NDV La 
Sota-infected groups at 7dpi and in IBV ck/CH/LDL/97I- and AIV H9- infected groups at 
12dpi may be responsible for the tracheal epithelium lesions at different stages of avian 
respiratory viruses infections. Further studies should focus on the exact roles of 
phosphorylated HSPB1 in the replication and pathogenesis of IBV, NDV and AIV H9. 

Viral infections are capable of induce oxidative stress associated with the activation of 
phagocytosis and release of reactive oxygen species which positively modulate immune 
activation, eliminate viral infection and immune-induced cellular injury [64]. In this study, 
antioxidant stress proteins including peroxiredoxin (PRDX)-1, 3, 4, 6, glutathione S- 
transferase (GST) and superoxide dismutase 1 (SOD1) (Table 2 and Supporting Information 
Table S2) were universally altered in all viruses-infected groups. Most of altered proteins 


were up-regulated, which may be related to relief of cellular oxidative stress induced by IBV, 
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NDV and AIV H9 infections. Moreover, more extensive changes were observed in AIV H9- 
infected group, suggesting that AIV H9 may activate more intense oxidative stress. Among 
these antioxidant stress proteins, PRDX1 and PRDX4 have been reported to protect RSV- 
induced oxidative damage to the nuclear cytoskeletal proteins in epithelial cells [65], while 
down-regulation of antioxidant enzymes (SOD, GST and PRDX) induced by RSV 
contributes to RSV pathogenesis [66]. In addition, PRDX1 can activate secretion of 
proinflammatory cytokines by interacting with TLR4 [67]. PRDX4 plays critical protective 
roles in inflammatory diseases [38]. Similarly, we considered that the up-regulation of 
PRDX1 and PRDX4 during IBV, NDV, and AIV H9 infections may activate host antiviral 
immunity to eliminate viral infection and provide protective cellular mechanisms to reduce 
oxidative stress and inflammatory damage. However, the potential roles of down-regulated 
oxidative stress-associated proteins in the pathogenesis of IBV, NDV, and AIV needed to be 
investigated further. 

In this study, trachea was chosen for global analysis of the host response because it is the 
primary target organ of avian respiratory viruses. The comparative analysis of changes in 
tracheal protein profile induced by IBV, NDV, and AIV H9 in vivo provides insight into 
mechanisms of global host response to different avian respiratory virus infections. However, 
the sites of IBV, NDV, and AIV infection are actually not completely consistent, IBV 
replicates in the epithelium of lower respiratory tract of chicken [68], NDV and AIV H9 
subtype induce lesions in the upper respiratory tract [18, 69]. The roles of different sites of 
respiratory tract in this study had not been effectively measured. Moreover, infections of IBV, 
NDV and AIV H9 subtype viruses can elicit inflammatory responses characterized by 
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recruitment of myeloid or lymphoid cells. Consequently, the cellular types of tracheal tissues 
harvested from mock and virus-infected chickens for proteomics analysis may differ due to 
immune cell infiltration. As a result, the changes revealed by proteomic analysis of some 
proteins such as LGALS1, PPIA and PRDX1 that could also be expressed in immune cells in 
the course of viral infection may be caused or partly caused by cellular heterogeneity 
occurred during the inflammatory response to viral infection (including infiltration by innate 
and adaptive immune cells). This needed to be further investigated. In addition, though a 
number of differentially expressed proteins that have potential roles in revealing the common 
or specific aspects of viral pathogenesis and virus-host interactions were identified, the 
precise roles of these proteins during viral infection had not proved in this study. 
Nevertheless, the similar and different protein expression profiles induced by different avian 
respiratory viruses after infection might be indicators of the essential response of infected 
target cells, provide insights into viral pathogenesis and host antiviral defense, and supply 
useful clues for the development of novel prevention or therapeutics strategies against these 


viruses. 
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Figure 1. The viral load was quantified using real-time RT-PCR. The viral loads in tracheal 
samples of three sacrificial birds from each group at each time point were determined. The 
average viral copy number per | ug RNA of each group was calculated through standard 
curve based on ten-fold dilution series of standard template with known concentration. The 
limit of detection of the real-time RT-PCR reaction was 107, 10° and 10° copies of template 
for IBV, NDV and AIV, respectively. (A) Kinetics of viral loads of IBV ck/CH/LDL/97I. (B) 
Kinetics of viral loads of IBV H120. (C) Kinetics of viral loads of NDV La Sota. (D) Kinetics 
of viral loads of and AIV H9. (E) The differences in viral kinetics among IBV 
ck/CH/LDL/97I-, H120-, NDV La Sota-, and AIV H9-infected groups. Data represent means 
of three biological replicates per group. Error bars indicate standard error of the mean, and 
dpi represent days post-inoculation. The symbol * indicate significant differences (* represent 


p <0.05 and ** represent p <0.01) between different groups. 
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Figure 2. The representative gel images showing 2D-DIGE maps labeled the relative 
locations of spots that displayed significant quantitative differential expression (p<0.05 and 
average ratio>1.2 or <-1.2) during different sampling times. (A) 3 dpi; (B) 7 dpi; (C) 12 dpi; 
and (D) 21 dpi. Representative 2D-DIGE images of IBV ck/CH/LDL/97I-, H120-, NDV La 
Sota-, and AIV H9-infected groups with their respective sampling points were represented as 
Supporting Information Figure S1. Fluor-labeled protein samples were separated on 24-cm 

pH 3-10 linear IPG strips, subjected to SDS-PAGE, and then gels were scanned on a 
Typhoon 9400 scanner, and image analysis was performed with Ettan™ DeCyder Software 

version v6.5. The identified protein spots were marked with arrows and labeled with the 
respective Match ID listed in Supporting Information Table S2. The green and red arrows and 
ellipses indicating HSPB1 and ANXA1 identified from separated spots. 
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Figure 3. Western blot confirmation of representative proteins in IBV ck/CH/LDL/97I-, 


H120-, NDV La Sota-, and AIV H9-infected chickens. (A) The immunoblot analysis of 


ANXAI, ANXA2, PRDX1, HSPB1 and phosphorylated HSPB1. 


(B) The relative 


densitometric intensity of ANXAI, ANXA2, PRDX1, HSPB1 and phosphorylated HSPB1 in 


immunoblot analysis. B-Tubulin was used as an internal control to normalize the quantitative 


data. 
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Figure 4. Gene Ontology annotation analyses of differentially expressed proteins according 
to their biological process categories. This classification was produced by using the 


GOSlimViewer tool of the Agbase database (http://www.agbase.msstate.edu/). (A) Biological 


process classification of all indentified proteins. (B) Comparison of biological process 
classification of differentially expressed proteins in IBV ck/CH/LDL/97]-, H120-, NDV La 


Sota- and AIV H9-infected groups and between IBV ck/CH/LDL/97I- and H120-infected 


Figure 4. 
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Figure 5. Transcript alteration of selected differentially expressed proteins in IBV 
ck/CH/LDL/97I-, H120-, NDV La Sota-, and AIV H9- infected chicken. (A) Comparison of 
mRNA levels of selected proteins differentially expressed at 3, 7, 12, and 21 dpi. Total RNA 

extracts were prepared from the tracheal tissues of sacrificial birds in control and all virus- 
infected groups at four time points (3, 7, 12, and 21 dpi), measured by real-time RT-PCR. 
Data represent means of three biological replicates per group. Error bars indicate standard 

error. Gene names refer to Supporting Information Table $6. Samples were normalized with 
the expression of the 28S ribosomal RNA gene in the control group. (B) Average ratios of 


selected proteins determined by 2D-DIGE using DeCyder Software. 
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Table 1. Summary of differentially expressed spots and proteins in IBV ck/CH/LDL/97I-, 
H120-, NDV La Sota- and AIV H9-infected groups and between IBV ck/CH/LDL/97I- and 
H120-infected groups 


IBV IBV NDV La AIV IBV 
ck/CH/LDL/971/ H120/Control — Sota/Control H9/Control ck/CH/LDL/971/ 
Control H120 
spots proteins spots proteins spots proteins spots proteins spots proteins 
Total” We 44 65 39 63 4] 86 48 63 38 
Up- 26 19 20 15 38 27 39 30 39 28 
regulated 
b) 
Down - 46 29 45 28 22 21 47 25 24 17 
pounce 


a) Total number of differentially expressed spots and protein. 
b) Number of up-regulated spots and proteins. 
c) Number of down-regulated spots and proteins. 
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Table 2. List of differentially expressed proteins commonly induced by IBV ck/CH/LDL/97I, H120, NDV La Sota, and AIV H9 infections 


Groups containing commonly Number of commonly expressed Gene name of commonly Biological process” 
expressed proteins” proteins expressed proteins 
IBV ck/CH/LDL/971/Control 12 akr carbohydrate metabolic process 10 
IBV H120/Control ENO1 generation of precursor metabolites and 
NDV La Sota/Control ANXAI energy 3 
AIV H9/Control protein complex assembly 3 
IBV ck/CH/LDL/971/ H120 ENO3 lipid metabolic process 1 
CKB transport 4 
GST response to stress 3 
HSPB1 cytoskeleton organization 3 
HBAA signal transduction 2 
cell proliferation 1 
HBAD catabolic process 3 
LMNA biosynthetic process 4 
PGAM1 small molecule metabolic process 4 
PRDX1 macromolecular complex assembly 3 
IBV ck/CH/LDL/971/Control 2 TPI1 reproduction | 
IBV H120/Control TTR carbohydrate metabolic process 3 
NDV La Sota/Control generation of precursor metabolites and 
AIV H9/Control energy | 
protein complex assembly 2 
lipid metabolic process 1 
catabolic process | 
biosynthetic process | 
cell differentiation | 
growth | 
small molecule metabolic process | 
anatomical structure development 1 
IBV ck/CH/LDL/971/Control 1 C21 - 
IBV H120/Control 
NDV La Sota/Control 
IBV ck/CH/LDL/97I/ H120 
IBV ck/CH/LDL/971/Control 6 LGALS1 immune system process 2 
NDV La Sota/Control DNA metabolic process 2 
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AIV H9/Control ERP29 transport | 
IBV ck/CH/LDL/971/ H120 EIFSA2 signal transduction 2 
APOBEC1 cell differentiation 2 
cellular nitrogen compound metabolic 
PRDX4 process 2 
PRDX3 anatomical structure development 2 
IBV ck/CH/LDL/971/Control 1 ARPC5 protein complex assembly 1 
IBV H120/Control cytoskeleton organization | 
AIV H9/Control macromolecular complex assembly | 
IBV ck/CH/LDL/971/ H120 
IBV H120/Control 1 EIFSA1 translation 5 
NDV La Sota/Control cellular protein modification process | 
AIV H9/Control cellular amino acid metabolic process 1 
IBV ck/CH/LDL/97I/ H120 transport 3 
biosynthetic process 5 
small molecule metabolic process | 
IBV ck/CH/LDL/971/Control 4 ADA immune system process | 
IBV H120/Control ACTA2 signal transduction | 
NDV La Sota/Control NDUFBIO aes process 1 
iosynthetic process 3 
PPIC cellular nitrogen compound metabolic 
process 5 
nucleobase-containing compound catabolic 
process | 
small molecule metabolic process 5 
IBV ck/CH/LDL/971/Control 2 NUMA cell morphogenesis 1 
IBV H120/Control SEPT2 cell cycle 1 
AIV H9/Control mitosis | 
signal transduction | 
cellular component assembly | 
anatomical structure formation involved in 
morphogenesis | 
anatomical structure development 1 
cell division | 
IBV ck/CH/LDL/971/Control 1 NDUFV2 - 
IBV H120/Control 
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IBV ck/CH/LDL/97I/ H120 


IBV ck/CH/LDL/971/Control 
AIV H9/Control 
IBV ck/CH/LDL/971/ H120 


IBV H120/Control 
NDV La Sota/Control 
AIV H9/Control 


IBV H120/Control 

NDV La Sota/Control 

IBV ck/CH/LDL/97I/ H120 
NDV La Sota/Control 

AIV H9/Control 

IBV ck/CH/LDL/971/ H120 


5 GAPDH 
LDHA 
PPIA 
PRDX6 
HADH 


3 ANXA2 
EPYC 
HEBP1 


1 PSMB2 


1 COL2A1 


immune system process | 
carbohydrate metabolic process 7 
generation of precursor metabolites and 
energy 2 

translation 1 

protein folding 1 

cellular protein modification process 2 
lipid metabolic process 2 

response to stress 1 

cytoskeleton organization | 

cell death 3 

catabolic process 3 

biosynthetic process 2 

small molecule metabolic process 3 
membrane organization | 

transport 2 

vesicle-mediated transport 1 

cell differentiation 2 

growth | 

anatomical structure development 4 
protein maturation | 

membrane organization | 

catabolic process 2 


cell adhesion 1 

cell death 1 

embryo development 3 

cell differentiation | 

extracellular matrix organization 1 
homeostatic process | 

anatomical structure formation involved in 
morphogenesis | 

anatomical structure development 18 
neurological system process 2 
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IBV ck/CH/LDL/971/Control 
IBV H120/Control 

IBV ck/CH/LDL/971/Control 
NDV La Sota/Control 


IBV H120/Control 
AIV H9/Control 


IBV H120/Control 

IBV ck/CH/LDL/97I/ H120 
NDV La Sota/Control 

IBV ck/CH/LDL/971/ H120 


1 THBS1 - 
3 aldolase C : 
C4 
OGN 
5) HAPLNI1 cell adhesion 1 
Nt Ss, 
cell proliferation 
Me) biosynthetic process 5 
cell differentiation | 
cellular nitrogen compound metabolic 
process 6 
small molecule metabolic process 6 
anatomical structure development 2 
1 ECI1 - 
1 GNB2L1 protein complex assembly 1| 


cellular protein modification process 2 
transport | 

cell cycle | 

signal transduction 4 

cell death 3 

catabolic process 2 

embryo development | 
vesicle-mediated transport 1 

cellular nitrogen compound metabolic 
process | 

nucleobase-containing compound catabolic 
process | 

growth | 

locomotion | 

small molecule metabolic process | 
anatomical structure development 1 
cell motility 1 

cell division 1 
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AIV H9/Control 3 COPE transport 3 
IBV ck/CH/LDL/97I/ H120 SOD1 catabolic process 1 
SRL vesicle-mediated transport 2 
cellular nitrogen compound metabolic 
process | 
nucleobase-containing compound catabolic 
process | 
small molecule metabolic process 1 
IBV ck/CH/LDL/971/Control 3 S100A11 DNA metabolic process 1 
ANXA8 cellular protein modification process 1 
signal transduction 3 
ree cell proliferation 1 
biosynthetic process 5 
cellular nitrogen compound metabolic 
process 4 
growth | 
chromosome organization 1 
IBV H120/Control 1 UQCRFS1 generation of precursor metabolites and 
energy | 
NDV La Sota/Control 2: rbf transport 6 
VDACI response to stress 1 
cell-cell signaling 2 
catabolic process | 
biosynthetic process | 
cellular nitrogen compound metabolic 
process 2 
nucleobase-containing compound catabolic 
process | 
small molecule metabolic process 2 
neurological system process 3 
transmembrane transport 3 
AIV H9/Control 4 RPSA - 


CTHRC1 
PSMA6 
TAGLN 
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a) Groups incubated different viruses (IBV ck/CH/LDL/97I, H120, NDV La Sota and AIV H9) containing commonly expressed proteins. 
b) Biological process of commonly expressed proteins in different viruses infected groups. This annotation was obtained using GOSlimViewer tool at Agbase 
database (http://www.agbase.msstate.edu/). — indicated no GO annotation was acquired. 
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Table 3. Proteins identified both in our study and previous reports in proteomes of the infected tissues or cells induced by AIV and IBV infections* 


*NDV was not included in this table because the proteome changes of cells or tissues induced by NDV have not yet been reported. 


IBV 
IBV IBV IBV IBV IBV (ck/CH/LDL/971 
AIV HS” HSN1” HINI” H9” (ck/CH/LDL/971)” ( 120)” (Beaudette)” — (Beaudette)” (H120)® Ps and P15)” 
chicken 
chicken human lung human AGS embryonic 
brain cells cells DF-1 cells vero cells trachea chicken trachea 
Tissues/cells chicken trachea (Ref. 8) (Ref. 10) (Ref. 11) chicken trachea chicken trachea (Ref. 12) (Ref. 13) (Ref. 14) (Ref. 15) 
Protein name 
(abbreviation) RPSA EIF5A2 RPSA RPSA ARPCS5 CKB ARPCS TPIL ARPCS5S 
LGALS1 EIFSA1 LGALS1 akr HBAA akr TTR akr 
CKB GST CKB ENOL HBAD ENOL CKB ENOI ENOI ENOI 
HSPB1 GAPDH HSPB1 LGALS1 NDUFV2 UQCRFS1 EPYC LGALS1 LGALS1 
LMNA HBAA LMNA/C ANXAIL NUMA ANXAIL HBAA ANXAIL ANXAIL ANXAI 
PPIA HBAD PPIA ERP29 OGN THBS1 HBAD ERP29 ERP29 
PRDX4 LDHA PRDX4 GST PPIA GST NUMA GST GST 
SEPT2 NUMA HSPB1 PRDX6 HSPB1 NME1 HSPB1 HSPB1 
ARPCS5S NMEI1 LMNA PGAM1 LMNA PGAM1 LMNA LMNA 
akr PRDX6 GAPDH ANXA8 ADA ECH GAPDH 
ENOL PGAM1 ADA HADH EIFSA1L HEBPI1 EIFSA1 EIFSA1 
ANXAI HEBP1 PHB NDUFB10 SEPT2 MYOZ1 PHB 
ANXA2 HADH SEPT2 PPIC ANXA2 NDUFB10 ANXA2 ANXA2 
ENO3 MYOZ1 LDHA PRDX4 ENO3 PPIC LDHA 
HAPLN1 PRDX1 PRDX1 APOBECI1 PRDX1 PSMB2 PRDX1 
TPIL PSMA6 PRDX3 THBS1 HAPLNI NDUFV2 PRDX3 
TTR APOBECI1 C21 TPIL C21 C21 
COPE CTHRC1 ACTA2 S100A11 ACTA2 ACTA2 
COL2A1 PRDX3 EIF5A2 TTR EIF5A2 
SOD1 SRL ALDOC C4 
ERP29 TAGLN ENO3 
EPYC TPIL 


a) Proteins identified in this study that induced by AIV H9 infection. 


b) Proteins identified both in our study and previous reports induced by AIV infection. 


c) Proteins identified in this study that induced by IBV ck/CH/LDL/97I and H120 infections. 


d) Proteins identified both in our study and previous reports induced by IBV ck/CH/LDL/97I and H120 infections. 
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